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INTRODUCTION 


A major  problem  facing  designers  and  manufacturers  of  high  speed 
turbine  engines  is  the  vibration  and  dynamic  loads  caused  by  rotor 
shaft  unbalance.  The  wider  range  of  rotational  speeds  and  the  use  of 
longer  shafts  in  multi-spool  engines  has  made  this  problem  more  acute. 
Moreover,  the  use  of  more  flexible  bearing  supports  in  controlling  the 
location  of  the  critical  speeds  allows  the  dynamic  forces  to  have  a 
greater  effect  upon  rotor  shaft  excursions. 

In  order  to  alleviate  this  problem  manufacturers  of  several  new 
turbine  engines  have  Incorporated  squeeze  film  dampers  in  the  rotor  shaft 
bearing  supports.  The  experience  of  these  manufacturers  indicates  that 
the  use  of  squeeze  film  dampers  significantly  reduces  engine  vibrations 
and  lowers  engine  rejection  rates  due  to  rotor  shaft  unbalance.  Experi- 
mental studies  by  Cooper  [1]  and  White  [2]  of  rotor  shafts  supported 
by  squeeze  film  dampers  have  shown  the  effect  of  the  squeeze  film  in 
decreasing  vibratory  loads  and  in  increasing  the  region  of  shaft  stability. 

It  ' us  also  been  shown  [3,4j that  external  damping  can  suppress  non-aynchronous 
whirl  at  super  critical  speeds.  A review  of  this  and  other  research  will 
be  made  later  in  this  section. 

The  squeeze  film  damper  configuration  is  relatively  simple.  A 
schematic  diagram  showing  the  operation  of  a squeeze  film  is  shown  in 
Figure  1.  A typical  installation  of  a squeeze  film  damper  is  shown  in 
Figure  2.  The  outer  race  of  the  rolling  element  bearing  is  fitted 
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loosely  in  the  bearing  support  structure.  The  radial  clearance  may 
range  from  .01  to  .0005  of  the  bearing  radius,  depending  upon  the  design. 
Radial  motion  of  the  bearing  outer  race  is  permitted  but  rotation  is 
prevented  using  some  type  of  mechanical  constraint.  The  annular  clearance. 
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Squeeze  Film 


between  the  bearing  outer  race  and  the  housing,  is  filled  with  oil. 

The  dynamic  forces  acting  on  the  bearing  will  result  in  an  orbital  motion 
of  the  bearing  outer  race.  This  motion  produces  corresponding  hydro- 
dynamic  forces  in  the  squeeze  film  which  oppose  the  shaft  motion  and 
rotor  dynamic  forces.  Thus  the  squeeze  film  acts  as  a complex  spring 
and  damper  system  which,  if  properly  designed,  will  significantly  reduce 
the  forces  transmitted  to  the  bearing  supports.  However,  it  has  also 
been  shown  [2,  5]  that  an  improperly  designed  squeeze  film  damper  can 
magnify  the  forces  transmitted  to  bearing  supports.  Thus  the  result 
is  a worse  situation  than  if  the  bearings  were  mounted  in  rigid  supports. 

In  order  to  correctly  design  a squeeze  film  damper  it  is  necessary 
to  predict  the  hydrodynamic  forces  in  the  squeeze  film  as  a function  of 
the  shaft  motion  and  pertinent  squeeze  film  parameters.  The  rotor 
dynamics  of  the  shaft  supported  by  squeeze  film  dampers  may  then  be 
predicted  and  the  squeeze  film  parameters  modified  to  obtain  the  desired 
damping  of  the  rotor  motion  and  forces.  These  parameters  include  the 
length-diameter  ratio,  the  eccentricity-clearance  ratio,  fluid  viscosity, 
and  supply  pressure.  The  effect  of  the  inlet  conditions,  cavitation  and 
end  seals  must  also  be  determined. 

In  order  to  predict  the  force  coefficients  of  the  squeeze  film  the 
pressure  distribution  throughout  the  film  must  be  determined.  Analytical 
and  experimental  studies  of  journal  bearings  have  resulted  in  several 
different  theories  for  predicting  the  pressure  distribution  in  journal 
bearings.  The  theoretical  model  of  squeeze  film  dampers  and  journal  bearings 
are  identical,  except  that  no  existing  journal  bearing  theory  takes  into 
account  the  boundary  conditions  produced  by  end  reals.  Since  end  seals 
reduce  or  eliminate  axial  flow,  it  may  be  surmised  that  the  "long  bearing" 
theory  would  fit  this  case  more  closely  than  the  "short  bearing"  theory. 
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The  theoretical  work  developed  for  journal  bearing^  may  be  applied  to 
squeeze  film  dampers  if  it  can  be  determined  that  the  boundary  conditions 
are  similar  or  produce  similar  results. 

It  was  the  purpose  of  this  study  to  make  an  experimental  investigation 
of  the  hydrodynamic  squeeze  film.  The  object,  in  designing  the  experi- 
mental rig,  was  to  duplicate,  as  closely  as  possible,  the  configuration 
found  in  turbine  engines,  while  keeping  the  measurements  as  simple  as 
possible. 

The  use  of  end  seals  in  the  test  rig  allowed  investigation  of  their 
effect  upon  the  pressure  distribution.  Many  squeeze  film  damper  designs 
and  some  journal  bearing  installations  contain  end  seals  to  minimize  end 
leakage,  yet  little  analytical  or  experimental  work  has  been  accomplished 
to  determine  their  effect  on  the  hydrodynamic  film  pressures. 

Initial  measurements  and  observations  showed  that  it  was  also 
desirable  to  study  the  effect  of  the  inlet  conditions.  The  inlet  design 
of  the  experimental  rig  and  some  squeeze  film  dampers  do  not  meet  the 
boundary  conditions  imposed  by  any  existing  theory.  Thus  the  effect 
of  varying  inlet  conditions  was  included  in  this  investigation. 

The  oil  supply  pressure  was  adjusted  so  as  to  maintain  a positive 
pressure  throughout  the  film.  Cavitation  was  thereby  prevented  and  a 
full  360  degree  film  was  developed.  Although  cavitation  does  occur  in 
practice  for  both  journal  bearings  and  squeeze  film  dampers,  the 
uncavitated  case  provides  the  basis  for  understanding  the  effects  of 
cavitation.  Later  research  will  cover  the  case  of  the  cavitated  film. 

The  experimental  rig  also  has  the  capacity  to  study  the  pressure 
distribution  for  the  case  of  the  eccentricity  being  a nonconstant  function 
of  time.  Little  or  no  research  has  been  carried  out  for  this  case, 
although  this  is  the  condition  found  in  any  journal  bearing  or  squeeze 
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film  damper  witn  a varying  load.  The  rotor  shaft  of  any  turboshaft  or 
turbo }et  engine  operated  in  a horizontal  position  will  whirl  about  an  equili- 
brium point  below  the  centerline  of  the  bearing  dampers,  due  to  the  effect  of 
gravity.  Inertia  forces  due  to  vehicle  maneuvers  can  deflect  the  orbit  center 
even  further.  Centering  springs  are  often  used  to  prevent  contact  with  the 
damper  bearing  wall,  but  these  are  made  as  soft  as  practical  to  keep  the  rigid- 
body  critical  speeds  low.  The  off  center  orbit,  even  if  circular,  produces 
both  non- zero  radial  velocities  and  non-constant  angular  velocities  of  the 
line  of  centers.  If  the  orbit  can  be  assumed  small,  it  is  common  practice 
to  transform  the  response  forces  from  the  rotating  coordinates  into  inertial 
(fixed)  coordinates  and  calculate  ‘'linearized”  stiffness  and  damping  coefficients 
for  a rotor  dynamics  analysis  of  the  motion  "in  the  small”. 

The  actual  size  of  the  whirl  orbit  depends  on  the  magnitude  of  rotor 
imbalance,  or  on  the  magnitude  of  any  destabilizing  self-excited  forces 
induced  by  the  motion.  Since  the  damper  force  response  is  nonlinear,  the 
orbit  size  is  a critical  factor  in  the  analysis.  The  type  of  analysis  to 
be  performed  (e.g.  linearized,  steady-state,  nonlinear,  transient)  should  be 
determined  by  the  primary  purpose  of  the  damper,  which  is  usually  one  of  the 
following: 

(a)  To  reduce  synchronous  response  to  rotor  unbalance  for  the 
expected  range  of  unbalance  in  production  engines,  thus 
reducing  rejections  due  to  excessive  vibration. 

(b)  To  increase  reliability  and  flight  safety  by  reducing  the 
response  to  sudden  and  abusive  unbalance,  such  as  caused 
by  the  losr;  of  a turbine  blade. 

(c)  To  suppress  nonsynchronous  whirl  and  self-excited  dynamic 
instabilities,  such  as  caused  by  internal  friction  in  the 
rotor-shaft  assembly. 
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Most  damper  and  rotor  dynamics  analysis  performed  to  guide  engine  design 
has  been  aimed  at  objective  (a).  This  objective  can  usually  be  met  satisfac- 
torily by  a design  based  on  a small  orbit  (linearized)  analysis,  as  described 
above,  or  on  a centered  circular  orbit  analysis  in  which  the  damping  a..d  stiff- 
ness coefficients  are  constant  in  the  rotating  coordinates. 

It  is  becoming  increasingly  clear  thet  more  effort  must  be  directed  toward 
objectives  (b)  and  (c) . By  definition,  a rotor  dynamics  analysis  to  meet  objective 
(b)  must  consider  large,  off  center,  whirl  orbits  with  non-zero  radial  veloc’ties. 
To  be  complete,  a stability  analysis  to  meet  objective  (c)  should  consider  both 
small  and  large  motions.  A transient  marching  solution  of  the  complete  equations 
by  digital  computer  is  the  only  method  capable  of  yielding  accurate  results  for 
some  of  these  cases,  but  the  results  can  be  no  more  accurate  than  the  force 
response  prediction  model  for  the  dampers. 

In  light  of  the  above,  a basic  objective  of  the  work  reported  here  was 
to  develop  an  experimental  capability  to  measure  damper  pressures  and  forces 
for  both  the  case  of  a constant  and  time  varying  eccentricity  (non-zero  radial 
velocity) , so  as  to  determine  the  validity  of  existing  theories  for  use  as  a 
prediction  model. 

As  will  he  described  in  a later  section,  the  resulting  experimental  rig 
produces  circular  orbits,  with  a variable  amount  of  offset  from  the  damper 
center.  It  is  felt  that  little  generality  is  lost  by  the  circular  constraint. 
Although  there  are  an  infinity  of  orbit  shapes  which  can  be  and  are  produced 
by  engine  rotors,  any  point  on  any  shape  orbit  can  be  uniquely  described  by  the 
instantaneous  eccentricity,  radial  velocity,  and  angular  velocity  of  the  line 
of  centers.  These  latter  quantities  are  all  reproducible  by  the  circular 
orbit  rig  with  the  proper  amount  of  offset. 
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ft  is  rupropriate  at  this  time*  to  briefly  survey  some  of  the  analytical 
and  experimental  research  that  hr.s  bean  accomplished  by  other  investigators 
related  co  this  study. 

Cooper  (11  was  one  of  the  first  investigators  to  study  the  effect  of 
the  bearing  supports  upon  the  rotor  dynamics.  In  an  experimental  study  he 
demonstrated  the  use  of  oil  films  in  controlling  shaft  vibrations. 

In  Cooper's  initial  tests  a stiff  shaft  was  supported  by  gas  bearings. 

The  effect  of  limiting  the  orbit  amplitude  using  mechanical  stops  was  studied 
first.  The  orbit  amplitude  was  then  limited  using  a In drodynamic  oil  film. 

It  was  found  that  limiting  the  journal  orbit  using  mechanical  stops 
caused  violent  pounding  of  the  stops  and  instability  of  the  shaft  motion. 
Introduction  of  the  hydrodynamic  oil  film  caused  smooth  running  of  the  shaft 
until  the  onset  of  oil-whirl  as  the  shaft  rotational  speed  reached  twice  the 
shaft  critical  speed. 

The  shaft  was  then  supported  using  rolling  element  bearings  mounted 
in  squeeze  film  dampers.  These  tests  resulted  in  smooth  running  of  the 
rotor  without  the  problem  of  oil-whirl.  These  tests  also  showed  the  exis- 
tence of  two  stable  orbits,  the  larger  being  characterized  as  a whirl  mede 
and  the  smaller  characterized  as  being  the  inverted  mode. 

White  (21,  extending  the  work  of  Cooper,  studied  the  dynamics  of  a 
simple  rigid  shaft  rxmnted  in  squeeze  film  dampers.  A theoretical  analysis, 
representing  the  squeeze  film  by  the  "short  bearing"  solution  of  Reynolds' 
equation,  considered  the  stability  cf  the  three  predicted  bearing  center  orbits. 

It  was  shown  that  the  intermediate  eccentricity  orbit  is  always  unstable 
and  that  the  whirl  (large)  and  inverted  (small)  orbits  are  always  stable. 

The  whirl  condition  was  found  to  be  undesirable  as  it  increased  the  forces 
transmitted  to  the  bearing  supports. 
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An  apparatus  designed  to  simulate  the  theoretical  model  was  constructed. 


Too  existence  of  c’.e  two  stable  equilibrium  orbits  was  observed.  It  was 
noted  that  there  were  speed  ranges  over  which  either  orbit  could  exist  for 
a pari icti\.r  out-of-halance.  It  was  clearly  shown  that  for  the  whirl  orbit, 
the  rotor  mass  center  was  more  distant  from  the  center  of  rotation,  thus 
the  squeeze  film  effectively  increased  the  unbalance. 

In  similar  research  Mohan  and  Hahn  {51  also  predicted  the  two  stable 
orbits  using  the  "short  bearing"  analysis  to  describe  the  squeeze  film.  The 
inverted  mode  was  characterized  by  a low  ...ansmissibility , less  than  one, 
while  the  whirl  mode  commonly  had  a transmissibility  greater  than  one.  A 
critical  shaft  unbalance  was  found,  above  which  only  the  whirl  mode  was  stable 
and  the  transmissibility  was  greater  than  one.  Design  graphs  were  presented 
to  aid  in  the  design  of  squeeze  film  dampers  over  a wide  range  of  operating 
comli  i ions. 

Vance  and  Lee  f3]  performed  a stability  analysis  of  high  speed  rotors 
with  internal  friction.  It  was  shown  that  the  threshold  speed  of  instability 
for  nonjyncbronous  whirl  induced  by  the  internal  friction  depends  on  the 
ratio  of  internal  friction  to  external  damping.  It  was  also  shown  that 
when  the  external  damping  exceeds  the  internal  friction,  operating  speeds 
up  to  eighty  percent  greater  than  the  critical  speed  may  be  safely  attained. 
Conversely*  when  the  internal  friction  is  many  times  larger  than  the  external 
damping  the  stability  speed  threshold  approaches  the  critical  speed.  It 
mav  he  inferred  that  correctly  designed  squeeze  film  bearing  dampers  can 
provide  the  external  damping  required  for  safe  operation  above  the  critical 
spued . 

Sweet  and  (Jenin  {A],  in  an  analytical  study,  demonstrated  that  the 
journal  hearing  instability  known  as  "oil  whip"  may  he  eliminated  hy  a 
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s. 


sufficient  amount  of  external  damping.  However,  the  magnitude  of  external 
damping  required  to  significantly  alter  the  undamped  results  is  much  too 
large  to  he  produced  by  normal  external  dissipative  forces.  It  was  then  shown 
that  a squeeze  film  damper  around  the  bearing  can  supply  the  necessary  external 
damping  to  eliminate  the  instability  due  to  oil  whip. 

The  above  studies  were,  for  the  most  part,  concerned  about  the  e'peet 
of  the  hydrodynamic  film  upon  the  dynamics  of  the  rigid  rotor.  The  nature 
of  the  hydrodynamic  film  was  not  studied.  The  theories  that  were  used  to 
predict  the  film  forces  have  not  yet  been  experimentally  verified. 

Thomsen  and  Anderson  [6],  in  an  experimental  study,  investigated  the 
hydrodynamic  squeeze  <r<1*n.  A vertical  shaft  was  used,  with  the  upper  rolling- 
clement  bearing  mounted  in  a centrally  preloaded  squeeze  film  damper.  The 
excitation  force  was  provided  by  deliberately  creating  an  unbalance  in  the 
shaft.  The  squeeze  film  damping  coefficient  was  obtained  by  measuring  the 
deflection  in  the  radial  supports  of  the  bearing  housing.  The  investigation 
was  primarily  intended  to  study  the  range  of  damping  available  from  the  snueeze 
film  b\  varying  radial  clearance  and  oil  viscosity. 

Potential  errors  are  incurred  in  determining  the  squeeze  film 
forces  by  measuring  the  deflection  of  the  housing  supports,  due  to  the 
effect  of  the  inertia  of  the  housing,  which  can  be  significant.  Also, 
deflection  of  the  housing  alters  the  squeeze  film  geometry,  introducing 
another  potential  source  of  error. 

Jones  [71  investigated  the  fundamental  hydrodynamics  of  the  squeeze 
film.  The  axial  and  circumferential  pressure  distribution  of  the  squeeze 
film  were  measured  about  a journal  constrained  to  move  in  a circular  orbit 
around  the  housing  center.  The  design  of  the  rig  was  such  as  to  approximate 
the  boundary  conditions  of  the  "short  bearing"  theory.  Oil  was  provided  to 
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the  squeeze  film  through  a central  circumferential  groove.  Cavitation  was 
allowed  in  the  divergent  part  ol  the  squeeze  film.  The  parameters  varied 
during  tne  tests  were  the  film  geometry,  journal  eccentricity  ratio,  orbital 
speed  and  the  supply  pressure. 

It  was  found  that  the  cavitated,  "short  bearing"  theory  provided  a 
reasonably  good  approximation  to  the  results  from  the  experiments.  Jones* 
study  gives  good  insight  into  the  squeeze  film  hydrodynamics  for  damper 
configurations  which  approximate  the  boundary  conditions  of  the  "short 
bearing"  theory. 

There  has  been  little  research  of  the  rotor  dynamics  or  squeeze  film 
hydrodynamics  for  configurations  approximating  the  "long  bearing*  theory. 

This  is  especially  true  for  the  case  of  a nor • constant  eccentricity  (i.e. 
non-zero  radia-  velocity). 

All  of  the  research  discussed  above  is  based  on  Reynolds’  theory  for 
laminar  flow.  Squeeze  film  dampers  may  operate  in  different  flow  regimes. 

Vohr  [Hi  investigated  the  fluid  flow  between  nonconcenti ic  rotating  cylinders 
(no  orbiting).  The  object  of  the  investigation  was  to  study  the-  nature  of 
Taylor  vortex  flow  and  turbulence  in  concentric  and  eccentric  annuli.  Cri- 
teria were  developed  as  a function  of  the  Reynolds'  number  and  Taylor  number 
for  the  onset  of  Taylor  vortices.  Although  the  results  may  not  be  quantita- 
tively applicable  to  squeeze  film  flow,  they  do  suggest  that  the  possibility 
of  flow  transition  is  a critical  question  to  be  considered. 


THEORY  OF  THE  HYDRODYNAMIC  ? QUEEZE  FILM 


The  equations  governing  the  hydrodynamics  of  the  squeeze  film 
were  developed  largely  fro',;  the.  theory  describing  journal  bearings. 
The  theoretical  model  and  sign  convention  are  shown  in  Figure  3. 

The  general  equation  for  the  hydrodynamic  film  between  two  flat 
surfaces  was  developed  by  Reynolds. 
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Implicit  in  the  above  equation  are  Reynolds'  assumptions  of  laminar  flow, 

incompressibility,  constant  viscosity,  and  that  the  Inertia  and  body 

0 

forces  are  negligible.  It  is  also  assumed  that  at  the  lubricated 
surfaces  the  fluid  has  the  velocity  of  the  surfaces  .aid  that  at  the 
film  extremities  the  pressure  depends  upon  the  external  conditions. 

It  should  be  noted  that  end  seals  do  not  meet  the  boundary  conditions 
imposed  by  Reynolds'  equation,  fhe  end  seals  prevent  the  pressure  at  the 
film  extremities  from  approaching  the  external  conditions.  The  end  seals 
also  constrain  the  velocities  at  the  extremities  to  equal  the  velocity 
of  the  seals. 

Reynolds'  equation  may  be  written  in  terms  of  Che  squeeze  film 
parameters , defined  in  Figure  3 and  the  List  of  Symbols . 


k ia  + n cos0^!q1  + rj2  fej t(1  + n cos0)3l^ 1 

121*; 

* {ursine  + n cos0] 


In  the  derivation  of  this  equation  it  is  assumed  that  the  effect  of  the 
curvature  is  negligible  and  certain  assumptions  are  made  concerning  the 
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relative  size  of  some  of  the  geometric  ratios. 

An  analytical  solution  to  Reynolds'  equation  has  never  been  obtained. 
Although  numerical  solutions  have  been  obtained , they  are  considered 
here  to  be  too  unwieldy  for  use  in  predicting  the  damper  forces  in  a 
useful  form  for  rotor  dynamics  analysis.  The  equations  can  be  consider- 
ably simplified  by  either  of  two  approximations,  both  of  which  Involve 
the  length  of  the  bearing  and  the  amount  of  fluid  flow  in  the  axial 
direction. 

If  the  fluid  film  is  assumed  to  be  infinitely  long  in  the  x^ 
direction  the  variation  of  the  film  properties  In  that  direction  will 
be  negligible.  Thus  the  second  term  on  the  left  hand  side  of  Reynolds' 
equation  may  be  neglected.  This  is  the  "long  bearing"  assumption  and 
is  the  basis  of  Sommerfeld’s  solution  for  journal  bearings. 

The  simplified  Reynolds'  equation  obtained  by  this  assumption  must 
then  be  integrated  twice  with  respect  to  9 to  obtain  the  pressure  distri- 
bution. This  integration  requires  the  choice  of  two  boundary  conditions. 

In  the  theory  developed  by  Sommerfeld  one  of  these  conditions  s'  acifles 
that  the  pressure  at  0-*0  be  a reference  pressure,  Pq.  For  the  case  of 
n radial  inlet  located  centrally  in  the  oil  film  the  validity  of  this 
assumption  should  be  questioned.  In  a journal  bearing  under  a constant 
load  the  inlet  position  would  not  vary  with  respect  to  0 and  the  assumption 
would  seem  to  be  valid.  In  a squeeze  film  damper,  however,  the  line  of 
centers  (the  line  passing  through  the  center  of  the  journal  and  the  center 
of  the  housing)  is  constantly  rotating.  Thus  the  position  of  the  inlet 
is  constantly  changing  with  respect  to  0.  The  effect  of  this  upon  the 
pressure  distribution  is  dependent  upon  the  ability  of  the  inlet  flow 
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to  dynamically  follow  the  changing  pressures  throughout  each  cycle.  A more 
elaborate  discussion  of  this  problem  is  contained  in  the  section  containing 
the  experimental  results. 

.rhe  result  of  the  second  integration  is  an  expression  for  the  circum- 
ferential pressure  distribution  about  the  journal.  Trumpler  [9],  using  this 
method  of  analysis,  obtains  the  following  expression  for  the  pressure 


distribution: 


bur.,  r / n x sin8(2+n  cos6)  , r 

v-Po  - (iVcc.-e/  + 1 


(1+n  ccs0)  ‘ 


1 1 x n-i 

(1+n  co80)  ^ (l+n)/J  n 


The  alternative  simplifying  assumption  concerning  Reynolds'  equation  is 
that  the  fluid  film  is  Infinitely  short  in  the  direction.  Thus  the  variation 
of  the  film  properties  in  that  direction  predominate  and  the  first  term  on 
the  left  hand  side  of  Reynold's  equation  may  be  neglected.  This  is  the  basis 
of  Dubois'  and  Ocvirk ' s p.01  solution  for  journal  bearings  and  is  known  as  the 
"short  bearing"  assumption. 

The  simplified  equation  obtained  using  this  assumption  is  then  integrated 
twice  with  respect  to  Xj  to  obtain  the  pressure  distribution.  Again,  this 
integration  requires  two  boundary  conditions.  These  specify  that  the  pressure 
at  each  end  of  the  film  is  equal  to  the  supply  pressure.  In  practice  the  oil 
is  usually  supplied  to  the  film  through  a central,  circumferential  groove  and 
the  ends  of  the  film  are  at  atmospheric  pressure.  Unless  the  lubricant  is 
supplied  at  near  atmospheric  pressure,  the  use  of  a central  supply  groove  does 
not  meet  the  specified  boundary  conditions.  It  is  evident  that  the  use  of  end 
seals  would  violate  these  boundary  conditions  to  an  even  greater  extent. 

The  result  of  this  integration  is  an  expression  for  the  pressure  distri- 
bution in  the  circumferential  and  axial  directions.  Jones  [7],  using  tL.s 


analysis,  obtains  (modified  for  the  sign  convention  used  in  this  study): 
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Neither  the  "long  bearing"  theory  nor  the  "short  bearing"  theory  is 
entirely  accurate  for  bearings  of  finite  length.  The  "short  bearing" 
assumption  is  usually  considered  valid  for  length/radius  ratios  less  than 
two  while  the  "long  bearing"  theory  is  considered  valid  for  1*.  ~!gth/radius 
ratios  greater  than  two.  Regard? ess  of  length,  a journal  bearing  or  squeeze 
film  damper  equipped  with  end  seals  may  approximate  the  "long  bearing"  theory 
in  that  flow  in  the  axial  direction  is  prevented. 

The  pressure  distribution  obtained  using  either  the  "short  bearing" 
assumption  or  the  "long  bearing"  assumption  may  be  integrated  to  obtain  the 
radial  and  tangential  force  components  of  the  squeeze  film. 

Using  the  "long  bearing"  theory  the  pressure  distribution  Is  integrated 
with  respect  to  6. 


pressure  distribution  is  integrated  with 


sinO  dx^d6 
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For  the  "short  bearing"  theory  the 
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These  integrations  require  a further  assumption  regarding  the 
circumferential  length  of  the  squeeze  film.  Both  solutions  to  Reynolds* 
equation  predict  significant  negative  pressures,  with  respect  to  the  supply 
pressure,  in  the  divergent  half  of  the  squeeze  film  (0  < 6 < n).  However, 
the  lubricant  film  cannot  support  negative  pressures,  as  vaporization  of 
the  lubricant  will  occur.  This  the  solutions  of  Reynolds’  equation  are  valid 
only  if  the  supply  pressure  is  sufficiently  high  to  prevent  cavitation.  For 
this  case  the  pressure  distribution  is  integrated  from  6*0  to  6*2tt  to  give 
the  force  components  for  the  full  ::ili>. 

If  the  supply  pressure  is  net  sufficiently  high  to  support  a full  film 
then  the  film  will  be  cavitated  ever  part  of  the  circumference.  The  pressure 
in  this  region  cannot  fall  below  the  vapor  pressure  of  the  lubricant  and  may 
be  assumed  to  equal  either  the  vapor  pressure  or  atmospheric  pressure.  For 
this  case  the  most  common  assumption  is  that  the  pressure  in  the  divergent 
part  of  the  film  is  equal  to  atmospheric  pressure  while  the  pressure  in  the 
convergent  part  of  the  film  (*r  < 0 < 2ir)  is  given  by  the  solution  to  Reynolds' 
equation.  Integration  using  these  limits  results  in  the  force  components  for 
the  cavitated  film. 

Expressions  for  the  force  components,  using  this  method  of  analysis, 
for  the  "long  bearing"  and  the  "short  bearing"  solutions  to  Reynolds' 
equation,  with  full  qnd  cavitated  films,  are  given  in  Table  1. 

Other  researchers  have  obtained  different  expressions  for  the  force 
components . 

Tondl[ll]  divides  the  squeeze  film  force  into  four  force  coefficients. 
The  force  components  are  presented  in  the  form: 

Fx  =H>uL  l(-2uc)  Fj(n)  + 2n  F*(r.)] 

Fy  - 6yL  r*  t(-2wc)  F^n)  + 2a  F,(n)] 
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TABLE  1.  Theoretical  Force  Components  for  the  "Long  Bearing"  and  "Short  Bearing"  Solutions  to  Reynolds  Equation 


The  preceding  form  ,*s  useful  as  it  separates  the  contributions  of  the  radial 
and  orbital  iot<nB  of  the  journal,  to  the  radial  and  tangential  force 
components.  In  his  analysis,  Tondl  calculates  the  force  components  directly 
without  first  determining  the  pressure  distribution.  However,  in  determining 
the  location  of  tbs  point  of  maximum  pressure,  Tondl  implicitly  makes  the 
assumption  of  zero  n.  This  would  invalidate  the  results  for  the  case  of  a 
non-zero  h. 

Hori  [12],  in  an  analysis  describing  journal  motion  during  oil  whip, 
also  formulates  force  components  which  would  be  applicable  to  a squeeze  film 
damper.  However,  certain  assumptions  are  made  concerning  surface  velocities 
which  are  true  only  for  n * 0. 

Lund  and  Saibel  [13]  divide  the  force  components  into  eight  coefficients 
in  cartesian  coordinates.  The  coefficients  arc  in  the  form  of  spring  and 
damping  constants: 
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Coefficients  in  the  above  form  allow  the  designer  to  simulate  the  dynamics 
of  the  rotor-bearing-damper  support  system  and  use  a numerical  marching 
solution  to  study  the  transient  response. 

However,  such  solutions  are  of  little  use  without  confidence  in  the 
theories  used  to  obtain  the  force  coefficients.  Current  literature  is  not 
consistent  in  the  development  of  equations  applicable  to  squeeze  film  dampers, 
especially  in  the  case  of  a non-zero  n.  The  basic  assumptions  upon  which 
Reynolds’  equation  is  based  should  be  questioned  for  many  squeeze  film  damper 


installations.  Vohr  [8]  developed  criteria  for  vortex  and  turbulent  flow 
between  nonconcentric  rotating  cylinders.  Applying  these  criteria  to  one 
squeeze  film  damper  installation  indicates  that  the  flow  is  turbulent,  thus 
"viking  an  analysis  based  on  Reynolds’  equation  invalid. 

The  high  rotational  speeds  and  the  low  viscosity  oil  used  in  many 
squeeze  film  damper  installations  make  Reynolds'  assumption  of  laminar  flow 
questionable.  A further  discussion  of  this  topic  is  carried  out  in  a later 
section. 
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EXPERIMENTAL  EQUIPMENT  AND  METHODS 


In  designing  the  experimental  rig  it  was  desired  to  duplicate  the 
configuration  of  squeeze  film  damper  designs  found  in  many  turbine 
engines , while  taking  advantage  of  the  theory  of  similitude  to  allow 
the  use  of  s larger  clearance  (to  minimize  the  effect  of  manufacturing 
tolerances)  and  to  allow  operation  at  lower  speeds. 

The  objective  of  the  design  was  to  provide  a journal  with  a known, 
constrained  motion  withip  an  annular  clearance  filled  with  oil.  In  this 
way  the  squeeze  film  could  be  studied  independently  from  the  dynamics  of  the 
rotor-bearing-eupport  system. 

It  was  decided  to  measure  the  pressure  distribution  in  the  oil  film 
Instead  of  measuring  the  resultant  film  forces.  The  pressure  measurements 
allow  a more  direct  check  of  the  existing  squeeze  film  theories  and 
the  film  forces  are  easily  calculated  from  the  pressure  distribution.  The 
film  forces  would  have  to  be  determined  by  measuring  the  housing  deflection 
and  it  would  be  difficult  to  determine  the  effect  of  the  housing  inertia 
on  these  measurements. 

Squeeze  Film  Apparatus 

A schematic  diagram  of  the  basic  rig  is  shown  in  Figure  4.  A photo- 
graph of  the  actual  hardware  is  shown  in  Figure  5. 

The  journal  is  mounted  on  the  eccentric  lobe  of  a stiff  shaft  using 
ball  bearings.  The  shaft  rotates  on  ball  bearings  mounted  in  rigid 
supports. 

The  outer  housing  is  supported  by  end  pistes  using  locational  fits 
in  order  to  allow  it  to  rotate  through  360  degrees.  The  housing  is  in- 
dexed in  10-degree  increments  and  a pin  is  inserted  in  order  to  keep  it 
from  rotating.  Holes  are  tapped  in  the  outer  cylinder  for  installation 
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URE  i.  Photograph  of  Expertise 


o!  pressure  transducers  and  a thermocouple.  These  are  in  direct  contact 
with  the  oil  film.  Rotation  of  the  outer  housing  allows  the  circumferential 
pressure  distribution  to  be  determined  using  one  pressure  transducer. 


The  journal  is  pinned  to  the  outer  housing  using  four  axial  pins. 

These  pins  allow  radial  and  orbital  motion  of  the  journal  but  prevent  rel- 
ative rotation  between  the  journal  and  housing. 

Oil  is  supplied  to  the  squeeze  film  through  an  orifice  located  cen- 
trally in  the  outer  housing,  140°  from  the  pressure  transducers  and  ther- 
mocouple. The  orifice  size  was  varied  on  several  tests  to  determine  the 
effect  of  differing  inlet  conditions  upon  the  pressure  distribution.  The 
oil  is  supplied  under  pressure  using  compressed  air  and  a pressure  regulator. 

Two  oils  of  differing  viscosity  were  used  during  the  experimentation. 
Viscosity  information  on  the  oils  is  given  in  Appendix  I. 

The  ends  of  the  journal  are  sealed  against  the  end  plates  using  *o’ 
rings  to  prevent  axial  flow  of  the  oil.  Leakage  past  the  'o’  rings  was 
minimal. 

The  eccentric  shaft  is  driven  by  an  electric  motor  using  a flexible 
coupling  to  prevent  axial  loads  on  the  journal.  The  rotational  speed  of 
the  motor  was  fixed  at  1765  revolutions  per  minute. 

The  eccentric  shaft,  hence  center  of  the  journal  orbit,  can  be  offset 
from  the  center  of  the  outer  housing  by  installing  shims  under  the  eccentric 
shaft  supports.  In  this  way  the  eccentricity  of  the  journal  orbit,  rela- 
tive to  the  center  of  the  housing,  varies  throughout  each  cycle  and  a known 
n is  obtained.  Separate  tests  were  carried  out  with  the  eccentric  shaft 
centered  and  offset  from  the  housing  center. 

Static  balance  of  the  eccentric  shaft  and  journal  was  accomplished  by 
the  addition  of  a counterweight  at  one  end  of  the  shaft. 

The  basic  dimensions  of  the  test  rig  pertinent  to  squeeze  film  calcu- 


lations are  given  in  Appendix  II. 
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I nst romentat ion 


Four  holes  are  tapped  in  the  outer  housing  for  installation  of  pres- 
sure transducers.  During  tne  tests  Viatran  Model  210  strain  gage  type 
pressure  transducers  were  used  with  ranges  of  0 to  100  psig  and  0 to  200 
psia.  The  natural  frequency  of  the  t 'ansducers  is  20  kilocycles/second  and 
the  response  time  is  500  microseconds  to  reach  90%  of  full  scale  output. 

The  pressure  transducer  was  calibrated  using  a dead  weight  tester  manufac- 
tured by  Amthor  Testing  Instrument  Company.  The  resulting  calibration 
curves  arc  given  in  Appendix  EL  The  transducer  is  connected  to  the  oil 
film  through  a short  cylindrical  tapping  in  the  outer  housing.  The  re- 
sulting dead  volume  was  assumed  to  have  a negligible  effect  on  the  dyna- 
mic pressure  readings.  To  obtain  the  circumferential  pressure  distribu- 
tion the  outer  housing  was  rotated  through  360  degrees  in  10  degree  incre- 
ments. 

The  pressure  transducer  was  powered  using  a ten  volt  battery  power 
supply.  The  output  from  the  transducer  was  displayed  on  a Type  561  B 
Tektronix  Oscilloscope  after  being  amplified  by  a Type  3A3  Dual-Trace 
Differential  Amplifier.  The  signal  was  displayed  versus  a time  sweep  gen- 
erated by  a Type  2B67  Time  Base. 

The  time  sweep  was  triggered  externally  using  the  amplified  signal 
from  a Type  MM0002  Bruel  and  Kjaer  magnetic  transducer.  The  signal  was 
amplified  using  a Bell  and  Howell  Model  1-165  DC  amplifier.  The  time 
sweep  was  triggered  each  time  the  lobe  of  the  eccentric  shaft  was  in  the 
vertical  up  position.  The  positioning  of  the  transducer  was  determined 
by  using  the  amplified  signal  to  trigger  a strobe  light.  The  instanta- 
neous position  of  the  shaft  at  the  time  of  triggering  could  then  be  deter- 
mined and  the  tranducer  moved  to  bring  the  lobe  to  the  wrticai  up  position. 
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A knowledge  of  the  oil  temperature  of  the  squeeze  film  was  desired 
in  order  to  make  an  accurate  determination  of  the  oil  viscosity.  This 


was  accomplished  by  mounting  a copper-constant.™  thermocouple  in  a pres- 
sure fitting  and  installing  it  in  one  of  the  holes  provided  for  the  pres- 
sure transducers.  Thus  the  thermocouple  is  in  direct  contact  with  the  oil 
film.  An  ice  poi-.i.t  was  used  for  the  reference  junction  and  the  output  was 
displayed  on  a Hewlett  Packard  Model  3490A  Multimeter.  The  thermocouple 
voltage  output  versus  temperature  agreed  with  standard  copper-cons tantine 
calibration  data. 

The  rotational  speed  of  the  eccentric  shaft  was  measured  by  inputing 
the  amplified  signal  from  the  magnetic  transducer  into  a Monsanto  Model 
103A  Counter  Timer. 

An  overall  view  of  the  test  rig  and  instrumentation  is  shown  in  Fig- 
ure 6. 

Method  of  Data  Analysis 

As  discussed  In  the  previous  section,  the  output  from  the  pressure  trans- 
ducer was  displayed  on  the  oscilloscope  screen  versus  a time  sweep.  The 
time  sweep  was  triggered  when  the  eccentric  lobe  of  the  shaft  was  in  the 
vertical  up  position.  The  distribution  displayed  on  the  screen  was  a vol- 
tage versus  time  curve.  In  order  to  get  meaningful  results  this  curve 
must  first  be  converted  to  pressure  versus  angular  rotation  of  the  eccen- 
tric shaft  and  ultimately  to  pressure  as  a function  of  the  squeeze  film 
position  angle  0. 

In  order  to  minimize  the  time  required  to  obtain  the  data  from  the 
oscilloscope  screen,  and  to  increase  the  accuracy  of  the  data  recording 
process,  a black  and  white  picture  of  the  distribution  on  the  oscilloscope 
was  taken  using  a 35  mm  camera.  After  pictures  were  taken  at  all  the 
desired  locations  of  the  outer  cylinder  the  film  was  developed. 
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holograph  of  Tost  Rig  anil  Instrumental- Ion 


The  negatives  were  then  cut  and  installed  in  plastic  slides.  These  slides 
were  put  into  a projector  and  displayed  on  a screen.  The  data  were  then 
recorded  from  the  projection  on  the  screen. 

The  time  sweep  of  the  oscilloscope  was  adjusted  so  that  each  full  cycle 
of  the  curve  covered  thirty-six  divisions  on  the  oscilloscope  screen.  Thus 
each  division  represented  ten  degrees  of  shaft  rotation. 

The  data  obtained  from  the  slide  were  the  output  from  the  pressure 
transducer,  in  millivolts,  versus  shaft  rotation,  with  readings  taken  at  ten 
degree  intervals.  The  voltages  werf  then  converted  to  pressures  using  the 
transducer  calibration  curves.  Thus  pressure  as  a function  of  shaft  rotation 
was  known. 

Experimental  data  were  recorded  for  two  cases:  (1)  the  eccentric  shaft 

centered  in  the  outer  housing  and  (2)  the  eccentric  shaft  offset  from  the 
center  of  the  outer  housing.  Henceforth  these  two  cases  will  be  referred  to 
as  the  centered  orbit  case  and  the  offset  orbit  case. 

For  the  centered  orbit  case  the  angle  of  shaft  rotation  may  be  related 
to  0 , where  0 is  the  angle  between  the  point  of  maximum  film  thickness  and 
the  pressure  transducer,  measured  in  the  direction  of  shaft  rotation.  Figure 
7 shows  the  geometric  relationship  needed  to  determine  0 from  the  angular 
position  of  the  eccentric  shaft  and  outer  housing.  The  pressure  p(0)  refers 
to  the  pressure  recorded  by  the  transducer  as  the  eccentric  shaft  rotates  which 
is  not  necessarily  the  same  as  p(0). 

For  the  centered  orbit  case  it  is  evident  that  the  pressure  distribution 
versus  *>  should  be  independent  of  the  position  of  the  pressure  transducer, 
provided  there  are  no  external  or  unknown  effects. 

For  the  offset  orbit  case  there  are  several  additional  variables  that 
affect  the  pressure  distribution.  For  thin  case  n,  n,  and  uc  are  no  longer 
constants,  but  are  functions  of  the  angular  position  of  the  shaft.  Data 
were  taken  at  each  of  the  36  positions  of  the  outer  housing  as  it  was  rotated 


Trigger 


0 * ir  + ifi0  - ip 


■<p  * Angle  of  eccentric  shaft 
from  vertical  up  position 
* 0)  t 

c 


FIGURE  7.  Geometric  Relationship  Between  0 and  the  .v. -gular  Position  of  the 
Eccentric  Shaft 


through  360  degrees  in  10  degree  increments.  In  order  to  find  the  pres- 
sure distribution  about  the  journal  in  one  particular  location,  the  pres- 
sure corresponding  to  that  journal  location  must  be  obtained  from  each 
of  the  thirty-six  photographs.  For  example,  if  it  was  desired  that  the 
pressure  distribution  about  the  journal  be  found  when  the  eccentric  Jobe 
had  rotated  180  degrees  from  the  vertical  up  position,  then  the  pressure 
corresponding  to  the  eighteenth  division  of  each  of  the  photographs  would 
be  recorded.  The  angular  relation  between  the  journal  position  and  the 
pressure  transducer  determines  the  value  of  6 for  each  of  the  readings. 

The  pressure  distributions  obtained  for  both  cases  were  compared  to 
the  existing  squeeze  film  theories.  The  pressure  distributions  were  then 
numerically  integrated  to  obtain  the  film  force  coefficients. 

A description  and  listing  of  the  computer  program's  used  in  the  data 
analysis  are  given  in  Appendix  IV. 


Experimental  Procedure 


In  order  to  facilitate  the  duplication  of  the  data  contained  in  this 
report  and  to  aid  future  experimental  studies  using  this  test  rig,  a step 
by  step  procedure  of  the  data  taking  process  will  be  outlined. 

1.  Turn  on  all  instrumentation  two  or  three  hours  before  taking 
data  in  order  that  they  may  become  thermally  stable. 

2.  Obtain  ice  and  water  in  a thermos  bottle  for  the  thermocouple 
reference  junction. 

3.  Check  the  oil  supply  in  the  holding  tank. 

4.  Rotate  the  outer  housing  to  bring  the  supply  orifice  to  the.  top. 

Remove  the  supply  hose  and  orifice.  Fill  the  squeeze  film  clear- 
ance with  oil  from  the  holding  tank  while  rotating  the  eccentric 
shaft  by  hand.  Reinstall  the  orifice  and  supply  line,  making  sure 
there  is  no  air  trapped  in  the  supply  line. 

5.  Aljust  the  pressure  regulator  of  the  compressed  air  supply  to  ob- 
tiiin  the  desired  supply  pressure. 

6.  Check  the  voltage  output  of  the  DC  power  source  and  adjust  the  rheostat 
to  obtain  10.0  volts. 

7.  Set  the  oscilloscope  time  sweep  to  trigger  on  an  external  negative 
signal  with  a positive  slope.  Run  the  experimental  rig  and  ad- 
just the  time  sweep  so  that  one  pressure  cycle  covers  36  divi- 
sions on  the  oscilloscope.  Adjust  the  sensitivity  of  the  oscil- 
loscope to  get  the  desired  amplitude  of  the  transducer  signal  on 
the  screen.  Record  the  voltage/division  setting. 

8.  After  the  rig  has  been  still  for  3 to  5 minutes  record  the  supply 
pressure  as  shown  by  the  oscilloscope. 
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Using  a 35  ma  camera  and  ASA  400  film,  the  shutter  speed  should 
be  set  at  1/15  second  and  the  f-stop  set  .?.t  between  2.8  and  4.0. 
Rotate  the  outer  cylinder  to  the  desired  location  and  Install  the 
locking  pin.  Run  the  rig  until  the  desired  oil  film  temperature 
is  reached  as  determined  by  the  output  of  the  thermocouple.  A 
picture  is  jn  taken  of  the  oscilloscope  screeen  and  a record  is 
made  of  th.  picture  number,  outer  cylinder  location  and  thermo- 
couple output. 

The  outer  cylinder  is  then  rotated  to  another  location  and  step 
10  is  repeated  making  sure  that  all  pictures  are  takrn  at  the 
same  oil  film  tetq>erature. 


RESULTS  AND  DISCUSSION 


Problems  with  the  experimental  apparatus  prevented  this  study 
from  being  as  extensive  as  had  heen  intended  from  the  onset.  Much 
of  the  effort  was  directed  tow*.  4 redesign  of  the  test  rig  to  prevent 
conditions  not  associated  with  squeeze  film  flow  from  altering  the 
pressure  distribution  about  the  orbiting  journal.  These  problems 
only  point  to  the  difficulty  in  trying  to  describe  squeeze  film 
bearing  damper  operation  under  actual  conditions,  where  the  damper 
is  interacting  with  the  dynamics  of  the  rotor-bearing-support  system. 
During  actual  operation  the  external  conditions  which  can  alter  the 
squeeze  film  are  more  prevalent  and  much  more  difficult  to  control. 

In  spite  of  the  problems  which  occurred  during  this  study, 
valuable  insight  was  gained  into  the  nature  of  the  hydrodynamics  of 
the  squeeze  film. 

Problems  with  the  Experimental  Rig 

For  the  centered  orbit  case,  the  journal  orbits  about  the  center 
of  the  outer  housing  with  a radius  equal  to  the  shaft  eccentricity. 

Due  to  symmetry,  the  pressure  distribution  with  respect  to  0 
should  be  independent  of  the  position  of  the  outer  housing  and  pressure 
transducer.  During  the  initial  running  of  the  rig  it  became  evident 
that  the  pressure  distribution  was  changing  as  the  outer  housing  was 
rotated.  Because  there  was  no  apparent  explanation  for  these  changes, 
a fairly  extensive  period  of  trouble  shooting  an-’  redesign  was  Initiated 
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to  determine  the  cause  of  the  changes  and  to  modify  the  test  rig 
accordingly.  The  following  is  a discussion  of  some  of  the  problems 
encountered  in  the  rig  and  a description  of  the  modifications  carried 
out  to  correct  them. 

During  each  pressure  cycle,  high  frequency  pressure  fluctuations 
of  large  amplitude  were  observed,  superimposed  on  the  pressure  distri- 
bution displayed  on  the  oscilloscope  screen.  These  pressure  fluctuations 
had  a shape  resembling  an  amplitude-time  plot  of  the  damped,  free 
vibration  of  a spring-mass  system.  A large  radial  pin  was  being  used 
to  prevent  relative  rotation  between  the  journal  and  the  outer  housing. 
Due  to  the  shape  of  the  pressure  fluctuations  it  was  believed  that  they 
were  caused  by  the  journal  and  shaft  vibrating  due  to  impact  with  the 
stiff  radial  pin.  The  radial  pin  was  therefore  removed  and  four,  small 
axial  pins  were  used  to  constrain  the  journal.  The  axial  pins  were 
designed  to  bend  under  loads,  thus  reducing  the  Impact  forces  on  the 
journal. 

A belt  drive  had  been  used  between  the  electric  motor  and  the 
eccentric  shaft.  Misalignment  between  the  belt  pulleys  resulted  in 
an  axial  load  on  the  eccentric  shaft  and  journal.  The  axial  positioning 
of  the  journal  was  controlled  by  the  compression  of  the  o'rings  used 
in  sealing  the  ends  of  the  journal.  An  axial  load  on  the  journal 
therefore  results  in  expansion  and  contraction  of  the  o'rings  which 
could  cause  a volumetric  change  within  the  squeeze  film,  resulting  in 
pressure  fluctuations.  The  belt  drive  was  replaced  by  a flexible 
coupling  designed  to  be  unable  to  transmit  axial  loads  between  the 
motor  and  eccentric  shaft. 


It  had  also  been  observed  that  the  inlet  conditions  had  a 
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significant  effect  upon  the  pressure  distribution  measured  by  the 
pressure  transducer.  The  pressure  distribution  could  be  altered  by 
varying  the  length  of  the  supply  hose  or  by  restricting  the  hose  in 
any  way.  In  order  to  get  meaningful,  repeatable  results  it  was  felt 
necessary  tc  make  the  pressure  distribution  independent  of  the  supply 
hose  and  tank,  except  for  the  magnitude  of  the  supply  pressure. 

Initially  a .110  inch  diameter  inlet  was  used  in  the  outer  housing  to 
supply  oil  to  the  squeeze  film  clearance.  An  orifice  was  designed  to 
be  installed  in  the  inlet.  It  was  desired  that  the  orifice  design  be 
such  that  the  supply  pressure  would  be  maintained  in  the  squeeze  film, 
but  that  the  pressure  distribution  would  not  be  altered  during  a pressure 
cycle  by  flow  through  the  inlet.  Experimentation  with  several  dif- 
ferent designs  resulted  in  an  orifice  with  a 1/64  inch  diameter  hole, 
two  inches  long. 

The  o' ring  grooves  at  each  end  of  the  journal  were  redesigned  so 
that  the  o' rings  were  a snug  fit  In  the  grooves.  This  was  done  to 
ensure  that  the  o' rings  could  not  move  in  the  grooves  causing  a vol- 
umetric change  in  the  squeeze  film  clearance. 

Measurements  were  made  of  various  parts  of  the  disassembled  rig 
to  determine  if  they  were  machined  to  drawing  specifications.  Special 
care  was  taken  in  measuring  the  eccentric  shaft  and  outer  housing 
supports,  as  they  control  the  relative  position  of  the  journal  within 
the  housiug.  The  rig  was  reassembled  and  the  eccentricity  of  the 
journal  was  measured  through  the  holes  provided  for  the  pressure 
transducers.  These  measurements  Indicated  that  the  outer  housing  and 
eccentric  shaft  were  concentric  within  .0005  inch. 

These  modifications  resulted  in  a large  improvement  in  the 
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observed  pressure  distributions.  The  large,  high  frequency  pressure 
fluctuations  during  each  cycle  were  no  longer  present  and  the  pressure 
distribution  appeared  to  be  independent  of  the  supply  hose  and  tank. 
Changes  still  occurred  in  the  pressure  distribution  as  the  outer 
cylinder  was  rotated  but  not  to  the  extent  seen  previously.  The  cause 
of  these  changes  is  net  known  but  nay  be  attributed  to  some  type  of 
volumetric  change  of  the  squeeze  film  clearance  or  a change  in  the 
squeeze  film  geometry  due  to  dynamic  loads  on  the  rolling  element 
bearings  and  supports. 

Centered  Orbit  Case 


As  discussed  in  the  previous  section,  some  changes  still  occurred 
in  the  pressure  distribution  as  the  outer  cylinder  and  pressure  trans- 
ducer were  rotated.  Due  to  these  changes,  data  is  presented  for  four 
perpendicular  positions  of  the  outer  cylinder.  The  angular  location 
of  the  pressure  transducer  for  these  positions  Is  0,  90,  180  and  270 
degrees  from  the  vertical  up  position,  in  the  direction  of  shaft  rotation. 

Data  were  recorded  using  both  SAE  10  and  SAE  50  oil.  In  order  to 
facilitate  comparison  with  the  existing  theory  the  results  were  put 
into  the  non-dimensional  form: 

p-Pp 

yr^wc/c2 

The  results  are  shown  in  comparison  to  the  "long  bearing"  and  "short 

4 

bearing"  theory  in  Figures  8,  9,  10  and  11  for  the  SAE  10  oil  and  in 
Figures  12,  13,  14  and  15  for  the  SAE  50  oil. 

The  general  shape  of  the  experimental  curves  follow  the  "long  bearing" 
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"Long  Bearing"  Theory 

2-5  in.  "Short  Bearing"  Theor 

4.813  in.  / \ Experimental  Data 


FIGURE  9 - Circumferential  Pressure  Distribution,  SAE  10  oil.  Transducer  Position  90 


"Long  Bearing"  Theory 


FIGURE  10  - Circumferential  Pressure  Distribution,— SAE  10  Oil,  Transducer  Position  180 


Experimental  Data 


FIGURE  11  - Circumferential  Pressure  Distribution,  SAE  10  Oil,  Transducer  Position  270 ' 


.062  in.  - - - - -Experimental  Data 


-12  L FIGURE  12*  Circumferential  Pressure  Distribution,  SAE  50  Oil,  Transducer  Position 


Long  Rearing"  Theory 


FIGURE  13  - Circumferential  Pressure  Distribution,  SAE  50  Oil,  Transducer  Position  90' 


Long  Bearing"  Theory 


-12  1“  FIGURE  14.  Circumferential  Pressure  Distribution  ,SAE  50  Oii,  Transducer  Position  180 


Long  Bearing"  Theory 


FIGURE  15.  Circumferential  Pressure  Distribution,  SAE  50  Oil,  Transducer  Position  270 


theory  with  fairly  good  agreement.  The  peak  positive  and  peak  negative  pres- 
sures occur  at  the  predicced  values  of  0.  However,  there  is  a large  variation 
in  the  amplitudes  of  the  pressures,  with  the  SAE  10  oil  producing  peak  to  peak 
pressures  53%  greater  than  predicted  while  the  SAE  50  oil  produced  peak  to 
peak  pressures  337.  less  than  predicted. 

Thus  the  measured  data  deviates  from  the  linear  dependence  on  viscosity 
as  predicted  by  Reynolds'  equate  n. 

Since  Reynolds’  equation  is  based  on  an  assumption  of  laminar  flow,  one 
possible  explanation  is  suggested  by  Vohr's  [8]  study  of  flow  between  non- 
concentric  rotating  cylinders.  In  this  investigation  the  nature  of  Taylor 
vortex  flow  and  turbulence  was  studied.  It  wa;3  shown  that  after  the  transi- 
tion from  laminar  flow  to  vortex  flow  the  shape  of  the  resulting  pressure 
distributions  is  nearly  unchanged  while  the  peak  to  peak  pressures  are  sig- 
nificantly larger  than  predicted  by  the  "long  bearing"  theory.  Criteria 
for  the  onset  of  Taylor  vortices  were  developed  using  a modified  Taylor  number 

/ 1/2  T « 41.2,  / 1/2  T « Vcp  / c - Re  / c 

cr.  cr.  — - 

M rj 

For  journal  bearings,  the  velocity  V was  defined  as  the  surface  velocity  of 
the  journal: 

V = rjUj 

In  squeeze  film  dampers  the  rotation  of  the  line  of  centers  (ojc)  produces 
the  viscous  pressure  distribution  defined  by  Reynolds’  equation  in  a rotating 
coordinate  system.  Since  the  Reynolds'  number  represents  the  dimensionless 
ratio  of  inertia  forces  to  viscous  shear  forces,  there  is  a question  as  to 
what  value  of  V is  appropriate  to  use  in  predicting  flow  transition  for  this 
case.  Although  no  experimental  data  exists  to  verify  transition  values  for 
the  case  of  a rotating  line  of  centers,  an  argument  can  be  made  for  using: 
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This  definition  for  V results  in  a value  of  80.2  for  the  SAE  10  oil 


and  a value  of  14.0  for  the  SAE  50  oil.  This  indicates  that  the  higher 
than  predicted  pressures  are  due  to  a transition  from  laminar  to  vortex 
flow. 

The  pressure  data  was  numerically  integrated  to  obtain  the  radial 
and  tangential  force  components,  and  also  the  damping  coefficient  B given 
by  Reference  6.  An  explanation  and  listing  of  the  program  used  to  accom- 
plish this  is  given  in  Appendix  4.  The  force  components  for  the  SAE  10  and 
the  SAE  50  oil,  the  damping  coefficient  B,  and  the  values  predicted  using 
the  "long  bearing"  theory  are  presented  in  Table  2.  As  expected  from  the 
pressure  distributions,  the  SAE  10  oil  had  values  approximately  35%  greater 
than  predicted,  while  the  SAE  50  oil  had  values  approximately  36%  less 
than  predicted  by  the  "long  bearing"  theory.  Thomsen  and  Anderson  [61,  in 
their  experimental  study,  found  that  the  damping  coefficient  was  significantly 
larger  than  predicted  by  the  theory  when  they  used  light  oil  (6  centipoise) 
and  large  clearances  in  their  test  apparatus.  It  is  under  these  conditions 
that  the  transition  from  laminar  to  vortex  flow  would  be  expected  to  occur. 
Using  the  criteria  given  above  and  the  data  given  by  Thomsen  and  Anderson 
a value  of  59.4  is  obtained  for  the  transition  parameter.  Thus  it  may  be 
surmised  that  the  larger  than  predicted  damping  coefficient  was  obtained 
due  to  a transition  from  laminar  to  vortex  flow.  Although  there  is  no 
conclusive  experimental  evidence, experience  in  industry  has  indicated  that 
there  may  not  be  the  linear  relationship  between  the  force  components  and 
viscosity  as  predicted  by  the  laminar  theory.  A possible  explanation 
is  that  the  dampers  may  be  operating  in  the  vortex  or  turbulent  flow  regime 
where  the  force  components  are  not  linearly  dependent  on  the  viscosity. 


EXPERIMENTAL  DATA 


TABLE  2:  FORCE  COMPONENTS  FOR  CENTERED  ORBIT  CASE 


r 

A calculation  of  the  Reynolds  number  for  the  damper  in  a recently  developed 
high  speed  turboshaft  engine  yields  a value  of  about  2500,  indicating 
turbulent  flow.  Engines  with  more  moderate  shaft  s peeds  may  have  dampers 
operating  in  either  the  vortex  or  laminar  flow  regime . Thus,  a determination 
of  the  proper  flow  regime  appears  to  be  an  important  prerequisite  to  an 
accurate  prediction  of  che  dynamic  force  response  for  engine  dampers. 

* 

deviation  of  the  results  for  the  SAF,  10  oil  from  the  "long  bearing" 
theory  may  be  explained  by  a transition  from  laminar  to  vortex  flow.  However, 
this  is  not  true  for  the  SAE  50  oil.  The  AE  50  oil  produced  peak  to  peak 
pressures  33%  less  than  predicted  by  the  "long  bearing"  theory  and  force 
t components  36%  less  than  predicted  by  the  "long  bearing"  theory.  The  devia- 

tion is  perhaps  not  surprising  due  to  the  differences  in  the  boundary  con- 
ditions of  the  theory  and  the  test  rig.  As  discussed  previously,  the  inlet 
of  the  test  rig  does  not  meet  the  boundary  condition  imposed  by  the  theory. 
However,  with  the  1/64  inch  orifice  it  is  believed  that  the  inlet  flow  cannot 

► 

dynamically  follow  the  changing  pressures  experienced  by  the  inlet  and  that 
the  pressure  at  0 = 0 stabilizes  to  some  constant  pressure. 

I The  use  of  end  seals  results  in  another  deviation  from  the  boundary 

i i 

conditions  for  Reynolds'  theory.  End  seals  prevent  axial  flow,  consistent 
with  the  "long  bearing"  theory,  but  impose  the  condition  of  zero  fluid 
velocity  at  the  ends.  This  is  a possible  explanation  for  the  lower  pres- 
sures produced  using  the  SAE  50  oil. 

Observations  with  pressure  transducers  located  away  from  the  central 
plane  of  the  outer  housing  indicated  that  peak  to  peak  pressures  did  not 
vary  in  the  axial  direction,  but  hat  there  was  approximately  a 10%  drop 
of  p0  from  the  center,  toward  the  end  of  the  journal.  These  results  are 
reasonably  consistent  with  the  "long  bearing"  theory  which  predicts  negligi- 
ble pressure  variation  in  the  axial  direction. 
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Effect  of  the  Inlet  Conditions 


In  the  initial  running  of  the  test  rig  it  was  observed  that  the  inlet 
conditions  altered  the  pressure  distribution  as  recorded  by  the  pressure 
transducer.  In  order  to  minimize  the  inlet  effects  the  1/64  inch  orifice 
was  installed  in  the  inlet  for  the  above  tests. 

It  was  then  decided  to  study  the  effect  of  a relatively  unrestricted 
inlet.  The  orifice  was  removed  leaving  the  .110  inch  supply  hole.  The 
tests  were  run  using  SAE  10  oil  and  data  were  recorded  for  two  perpendicular 
positions  of  the  pressure  transducer.  The  results  are  shown  in  comparison 
to  the  "long  bearing"  theory  in  Figures  16  and  17.  It  should  be  emphasized 
that  these  data  do  not  represent  the  instantaneous  pressure  distribution 
about  the  journal,  but  are  the  pressures  recorded  by  the  pressure  transducer 
as  the  eccentric  shaft  rotates. 

With  the  SAE  10  oil  and  the  .110  inch  supply  hole  the  flow  through  the 
inlet  is  relatively  unrestricted.  Thus  the  pressure  of  the  squeeze  film  at 
the  location  of  the  inlet  remains  at  the  supply  pressure.  As  the  eccentric 
shaft  rotates,  the  position  of  the  inlet  with  respect  to  0 is  constantly 
changing.  If  it  is  assumed  that  the  Sommerfeld  distribution  is  still  valid 
then  the  value  of  p^  must  be  constantly  changing  in  order  that  p(0)  at  the 
inlet  equals  the  supply  pressure.  Thus  the  instantaneous  pressure  distribu- 
tion about  the  journal  would  be  given  by  the  Somraerfeid  theory  with  the 

value  of  p such  that  the  pressure  at  the  inlet  equals  the  supply  pressure, 
o 

Figure  18  shows  the  effect  of  this  condition. 

Under  the  above  assumptions  a simplified  theory  was  developed  to  predict 
the  pressure  distribution  as  recorded  by  the  pressure  transducer.  The  supply 
pressure  for  both  the  restricted  and  the  unrestricted  case  was  45  psig. 

Using  the  data  from  the  restricted  case  the  quantity  [45  - p(9^.)]  was  calculated 
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FIGURE  16.  Circumferential  Pressure  Distribution,  Unrestricted  Inlet,  Transducer  Position 


/ 
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Circumferential  Pressure  Distribution,  Unrestricted  Inlet,  Transducer  Position  ISO 


at  ten  degree  increments  of  0^,  where  0^  is  the  inlet  position  (see  Figure 
19).  The  quantity  (45  - p(0j)l  represents  the  amount  which  p(0)  must  be 
adjusted  to  predict  the  pressure  distribution  for  an  unrestricted  inlet. 

p(0)  unrestricted  * p(0)  restricted  + (45  - p(0^)  restricted] 

The  resulting  curves  are  shown  with  the  experimental  results  in  Figures  20 
and  21. 

As  can  be  seen  from  the  figures  there  is  good  agreement  between  the 
shapes  of  the  curves  resulting  from  the  simplified  theory  and  the  experimental 
results.  It  is  evident,  however,  that  the  inlet  is  not  unrestricted  and 
that  the  inertia  of  the  fluid  is  not  negligible.  Tt  is  most  probably  these 
effects  which  cause  the  phase  and  amplitude  differences  between  the  theory 
and  the  experimental  results. 

Tliis  would  indicate  that  the  simplified  theory,  and  the  assumptions 
upon  which  it  is  based  are  essentially  correct.  A radial  inlet  in  an  actual 
squeeze  film  damper  design  would  probably  fall  between  the  restricted  and 
unrestricted  inlets  studied  in  this  investigation.  However,  the  cyclic 
variation  of  the  pressure  distribution,  as  shown  in  Figure  18,  would  be 
present  to  some  extent,  with  a resulting  variation  of  the  force  components. 

If  the  supply  pressure  is  all  compared  to  the  peak  positive  and  peak  negative 
pressures,  as  is  true  for  actual  squeeze  film  dampers,  then  the  extent  of  the 
cavitated  film  region  will  vary  throughout  each  cycle.  This  may  result  in  the 
film  being  completely  cavitated  over  portions  of  each  cycle,  resulting  in  a 
cyclic  loss  of  the  load  carrying  ability  of  the  squeeze  film. 

Offset  Orbit  Case 

The  eccentric  shaft  was  offset  from  the  center  of  the  outer  housing 
by  installing  shims  under  the  eccentric  shaft  supports.  In  this  way  the 


52 


Simpliied  Theory 


FIGURE.  20.  Circumferential  Pressure  Distribution,  Unrestricted  Inlet,  Transducer  Position 


Simplified  Theory 


FIGURE  21.  Circumferential  Pressure  Distribution, Unrestricted  Inlet,  Transducer  Position  180 


journal  was  made  to  orbit  about  a point  above  the  center  of  the  outer 


housing  and  a known  n,  n,  and  u>c  were  obtained  as  functions  of  the  shaft 
position.  SAE  50  oil  was  used  for  these  tests. 

The  pressure  distribution  around  the  journal  was  found  by  recording 
the  pressure,  corresponding  to  that  particular  journal  location,  from  each 
of  the  photographs  for  the  36  positions  of  outer  housing.  As  discussed  in 
the  previous  section,  for  the  centered  orbit  case,  changes  were  occurring 
in  the  pressure  distribution  as  the  outer  housing  was  rotated.  The  effect 
of  these  changes  on  the  accuracy  of  the  data-taking  process  for  the  offset 
orbit  case  cannot  be  determined,  but  may  be  significant. 

As  the  outer  housing  was  rotated,  the  pressure  distribution,  as  exper- 
ienced by  the  inlet,  was  constantly  changing.  Due  to  this,  the  value  of  pQ 
for  each  different  position  of  the  outer  housing,  corresponding  to  each 
photograph,  was  different.  To  correct  for  this,  a specific  value  of  pQ  was 
chosen  and  each  pressure  distribution  for  the  36  positions  of  the  outer 
cylinder,  was  adjusted  so  that  each  distribution  had  the  same  value  of  pq. 

Data  are  presented  for  four  perpendicular  positions  of  the  eccentric 
shaft.  These  positions  are  0,  90,  ISO  and  270  degrees  from  the  vertical  up 
position,  in  the  direction  of  shaft  rotation.  The  experimental  results  are 
shown  in  Figures  22,  23,  24  and  25  in  comparison  to  the  "short  bearing"  and 
"long  bearing"  solutions  to  Reynolds’  equation. 

As  can  be  seen  from  the  rigures,  there  is  considerable  difference  between 
the  experimental  results  and  the  theory.  In  Figures  22  and  24  the  eccentric 
shaft  Is  in  the  vertical  up  and  vertical  down  position  where  n is  zero.  Thus, 
in  looking  back  to  the  results  of  the  centered  orbit  case,  where  n is  also 
zero,  it  would  be  expected  that  the  experimental  results  should  follow  the 
theory  wiih  reasonable  agreement,  at  least  to  the  general  shape  of  the  curves. 
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Short  Bearing"  Theory 


Long  Bearing”  Theory 


Short  Bearing"  Theory 


FIGURE  25.  Circumferential  Pressure  Distribution,  Offset  Eccentric  Shaft  .osition  270 


This  is  true  only  for  certain  portions  of  the  curve  as  seen  . Figure  24. 

This  would  Indicate  that  either  the  measurements  are  inaccurate  or  that  the 
fluid  flow  characteristics  are  altered  because  the  journal  has  a non-zero 
h and  <i)c  over  all  other  portions  of  its  orbit. 

The  pressure  distribution  about  the  journal  when  located  90  and  270 
degrees  from  the  vertical  up  position  are  shown  in  Figures  23  and  25.  For  these 
two  positions  n and  Wc  are  identical  and  h is  of  the  same  magnitude  but  differs 
in  sign.  In  both  cases  the  experimental  results  do  not  follow  the  theory.  It 
should  be  noted  that  for  both  cases  a reasonably  smooth  curve  can  be  fitted 
through  the  experimental  results. 

The  pressure  data  was  numerically  integrated  to  obtain  the  radial  and 
tangential  force  components  and  the  damping  coefficient  B.  The  force  components, 
damping  coefficient  B,  and  the  values  predicted  using  the  "long  bearing"  theory 
are  presented  in  Table  3.  As  would  be  expected  there  is  considerable  disagreement 
between  the  force  components  derived  experimentally  and  those  predicted  by  the 
"long  bearing"  theory. 

As  discussed  previously,  the  analytical  treatment  in  current  literature 
of  the  case  of  a non-zero  ft  is  not  consistent  for  either  squeeze  film  dampers 
or  journal  bearings.  There  is  also  no  experimental  verification  of  any  theory 
for  the  case  of  a non-zero  h.  Thus,  the  results  of  this  study  should  not  be 
discourted  simply  because  they  do  not  agree  with  the  existing  theories.  The 
large  disagreement  between  these  results  and  the  "long  bearing"  theory  may 
indicate  that  the  "long  bearing"  theory  does  not  accurately  describe  the  squeeze 
film  for  the  case  of  a non-zero  ft. 

The  transition  from  laminar  to  vortex  or  turbulent  flow  has  not  been 
investigated  for  journal  bearings  or  squeeze  film  dampers  for  the  case  of  a 
non-zero  n.  It  is  possible  that  some  combination  of  ft,  wc  and  n over  portions 
of  the  journal  orbit,  for  the  offset  orbit  case,  result  in  the  formation  of 
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FORCE  COMPONENTS  FOR  OFFSET  ORBIT  CASE 


vortex  or  turbulent  flow.  This  would  account  for  the  discrepancies  between 
the  "long  bearing"  theory  and  the  experimental  results. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 

The  following  conclusions  were  drawn  from  this  study: 

1.  For  a damper  with  end  seals  preventing  axial  flow,  the  "long  bearing" 
theory  is  reasonably  accurate  in  predicting  the  shape  of  the  pressure 
distribution  for  the  centered  orbit  case.  This  holds  true  for  the 
vortex  flow  regime  as  well  as  for  laminar  flow.  However,  the  peak  to 
peak  magnitude  of  the  pressure  distribution  is  not  accurately  predicted 
by  this  theory,  even  in  the  laminar  flow  regime. 

2.  The  SAE  10  oil  produced  pressures  and  resulting  force  components 
significantly  greater  than  predicted  by  the  "long  bearing"  theory. 
Calculation  of  the  transition  parameter  presented  by  Vohr  [8]  indicates 
that  the  damper  was  operating  in  the  vortex  flow  regime,  resulting  in 
the  greater  magnitude  of  the  peak  to  peak  pressures.  Subsequent  calcula- 
tion of  this  parameter  for  actual  squeeze  film  dampers  indicates  that 
many  dampers  may  be  operating  in  the  vortex  or  turbulent  flow  regime. 

3.  The  SAE  50  oil  produced  pressures  and  resulting  force  components  signif- 
icantly less  than  preuicted  by  the  "long  bearing"  theory.  Calculation 
of  the  transition  parameter  indicates  that  the  damper  was  operating  in 
the  laminar  flow  regime.  A possible  explanation  for  the  deviation  of  the 
experimental  results  from  the  "long  bearing"  theory  is  the  differing 
boundary  conditions  of  the  test  rig  from  the  theory  due  to  the  use  of 
end  seals. 

A.  A radial  inlet  with  no  circumferential  grooving  results  in  significant 
deviation  from  the  "long  bearing"  theory.  An  unrestricted  radial  inlet 
results  in  a cyclic  variation  of  the  force  components  of  the  squeeze  film 
and  can  result  in  the  loss  of  the  load  carrying  ability  of  the  film  over 
portions  of  each  cycle. 
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Efficient  end  seals  (minimal  end  leakage)  result  in  a negligible  pressure 
variation  in  the  axial  direction,  consistent  with  the  "long  bearing"  theory. 
Although  the  experimental  results  for  the  offset  orbit  case  are  not 
sufficiently  accurate  or  complete  to  form  a reliable  prediction  model,  it 
can  be  stated  with  reasonable  confidence  that  the  "long  bearing"  theory 
does  not  accurately  predict  the  pressure  distribution  for  this  case  (cyclic 
variation  of  n,  n,  and  toc) . 


Recommendations  for  Future  Research 


It  Is  recommended  that  several  basic  modifications  * carried  out  on  the 
test  apparatus  to  extend  the  investigation  begun  in  this  study.  The  following 
is  a list  of  possible  modifications: 

1.  Machine  a central  circumferential  groove  in  the  outer  housing  to  be  used 
as  the  oil  inlet  to  the  squeeze  film.  This  will  eliminate  the  problems 
inherent  with  a radial  inlet  and  be  more  consistent  with  designs  currently 
being  used  by  several  engine  manufacturers. 

2.  Machine  the  outer  housing  for  additional  pressure  transducers  in  the 
axial  direction  to  determine  the  axial  pressure  distribution,  which  would 
be  more  significant  with  a central  circumferential  inlet.  The  addition 
of  several  thermocouples  can  determine  if  there  is  a significant  axial 

or  circumferential  temperature  gradient. 

3.  Replace  the  o'ring  seals,  used  at  the  journal  ends,  with  metal  piston 
ring  type  seals.  This  will  eliminate  the  possible  volumetric  change 
within  the  squeeze  film  clearance  cause  by  expansion  and  contraction  of 
the  flexible  o* rings. 

With  the  above  modifications,  the  centered  orbit  case  could  be  studied 
more  thoroughly.  Perhaps  one  of  the  most  interesting  areas  to  be  investigated 
is  the  relationship  between  lubricant  viscosity  and  the  magnitude  of  the 
resulting  squeeze  film  pressures.  This  study  could  also  lead  to  an  accurate 
determination  of  the  transition  from  laminar  to  vortex  or  turbulent  flow  in 
the  case  of  a rotating  line  of  centers  (orbiting  damper  journal). 

A more  accurate  method  of  determining  the  angular  position  of  the  eccentric 
shaft  is  needed  to  study  the  offset  orbit  case.  Data  for  this  case  consists 
of  discrete  data  points,  for  a particular  journal  location,  taken  from  each  of 
the  pressure  distributions  obtained  from  the  36  positions  of  the  outer  housing. 
This  requires  greater  accuracy  than  was  available  with  the  method  used  in  this 
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study.  Perhaps  high  speed  strip  chart  recorders  could  be  used  to  more 
accurately  record  this  data. 
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APPENDIX  I 


0 L VISCOSITY  INFORMATION 

Two  oils  of  differing  viscosity  were  used  during  the  experimen- 
tation. Viscosity  information  on  the  oils  was  obtained  from  the  manu- 
facturers. The  viscosity  and  related  information  for  the  two  oils 
is  given  below: 

SAE  10  Oil 

The  SAE  10  oi 1 was  produced  by  the  Ashland  Oil  Company  (Valvoline) . 
Viscosity  information  on  this  oil  was  obtained  from  the  Product  Appli- 
cations Department: 

Temperature  Viscosity 

100°F  186.7  SUS 

210°F  45.91SUS 

These  values  were  plotted  on  an  ASTK  (American  Society  for  Testing 
Materials)  Standard  Viscosity-Temperature  Chart.  When  the  Saybolt 
Universal  Viscosity-Temperature  relationship  is  plotted  on  the  ASTM 
chart,  the  resulting  curve  is  a straight  line.  Consequently  if  the 
viscosity  temperature  relationship  is  known  at  two  points,  i.e.  100°F 
and  210°F,  then  all  points  will  fall  on  the  line  connecting  these 
points.  Thus  the  viscosity  at  any  temperature  can  be  read  from  this 
line.  All  data  using  the  SAE  10  oil  vere  recorded  at  an  oil  temperature 
of  103°F.  The  viscosity  corresponding  to  this  temperature  is  170  SUS. 

Using  standard  unit  conversions  the  viscosity  was  obtained  in  reyns 
u * (4.6) (10  **) reyns  (lb-sec/in^) 
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SAE  50  Oil 


The  SAE  50  oil  was  produced  by  the  Standard  Oil  Company.  Viscosity 

information  was  obtained  from  them  on  request: 

Temperature  Viscosity 

100°F  1650  SUS 

210°F  90  SUS 

These  values  were  also  verified  by  another  researcher  at  the  University 
of  Florida. 

Again  these  values  were  plotted  on  an  ASTM  Standard  Viscosity- 
Temperature  Chart.  All  the  data  using  the  SAE  50  oil  were  recorded  at 
an  oil  temperature  of  115°F.  The  viscosity  corresponding  to  this 
* temperature  is  950  SUS. 

Using  standard  unit  conversions  the  viscosity  was  obtained  in  reyns: 
| y * (2.74)10  ^reyns  (lb-sec/in 
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DIMENSIONS  OF  SQUEEZE  FILM  APPARATUS 


ui 

(diameter  of  journal) 

5.0  in. 

rJ 

(radius  of  journal) 

2.5  in. 

L 

(length  of  journal) 

4.813  in. 

Db 

(inside  diameter  of  outer  housing) 

5.125  in. 

c 

(clearance  (Db-D^)/2) 

.062  in. 

Centered  Orbit  Case 

n 

(e/c) 

.50 

• 

n 

(e/c) 

0.0 

b) 

c 

(angular  speed  of  eccentric  shaft) 

184.9  rad/sec 

Offset  Orbit  Case 

A .015  inch  shin  was  placed  under  eccentric  shaft  supports.  Thus  the 
journal  was  made  to  orbit  about  a point  .015  in.  above  the  center  of 
the  outer  housing.  The  parameters  listed  below  varied  continuously 
between  the  given  values,  as  a function  of  the  eccentric  shaft  position. 


n 

(maximum) 

.742 

n 

(minimum) 

.258 

0) 

c 

(maximum) 

358.3  rad/sec 

jj 

c 

(minimum) 

124.6  rad/sec 

• 

n 

(maximum) 

44.7  1/sec 

. 

n 

(minimum) 

-44.7  1/sec 
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PRESSURE  TRANSDUCER  CALIBRATION  INFORMATION 

Two  Viatran  Model  210  strain  gage  type  pressure  transducers  were 
used  during  this  study.  Calibration  data  were  supplied  by  the  manu- 
facturer but  proved  to  be  inaccurate.  The  transducers  we-*e  therefore 
recalibrated  using  a dead  weight  tester  manufactured  by  Amthor  Testing 
Instrument  Company.  The  resulting  calibration  curves  are  shown  In 
Figures  26  and  27. 

The  calibration  data  were  also  used  to  derive  equations  to  be 
used  in  the  computer  programs  for  data  analysis.  The  resulting 
equations  are  given  below: 

Model  PTB710G-100  (Serial  No.  142571) 

Pressure  (psig)  « 5.5033  x voltage  output  (mv)  + .92 

Model  PTB210A-200  (Serial  No.  288074) 

Pressure  (psig)  ■ 8.146  x voltage  output  (mv)  - 11.349 


i 


100  r VIATRAN  PRESSURE  TRANSDUCER 


FIGURE  26.  Transducer  Calibration  Curvt  , 0 - 100  PSIG 


VIATRAN  PRESSURE  TRANSDUCER 
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DESCRIPTIONS  OF  COMPUTER  PROGRAMS  USED  IN  DATA  ANALYSIS 


Several  computer  programs  were  written  in  order  to  analyze  the 
data  obtained  from  this  study.  The  following  is  a description  of 
each  of  the  programs.  Included  is  a listing  of  each  program. 

1.  Calculation  of  the  Pressure  Distribution  Using  the  "Long  Bearing" 
Theory. 

This  program  calculates  the  circumferential  pressure  distribution 
about  a moving  journal  using  the  "long  bearing"  theory.  The 
journal  position,  radial  velocity,  tangential  velocity,  oil 
viscosity  and  other  pertinent  squeeze  film  parameters  are  input. 

The  output  is  the  pressure  distribution  about  the  journal  as  a 
function  of  0 , given  in  ten  uegree  increments. 

2.  Calculation  of  the  Pressure  Distribution  Using  the  "Short  Bearing" 
Theory. 

This  program  calculates  the  circumferential  pressure  distribution 
about  the.  center  of  the  journal  using  the  "short  bearing"  theory. 

The  journal  position,  radial  and  tangential  velocity,  oil  viscosity, 
and  other  pertinent  squeeze  film  parameters  are  input.  The  output 
is  the  pressure  distribution  about  the  center  of  the  journal  given 
in  ten  degree  increments  of  0 . 

3.  Calculation  of  the  Pressure  Distribution  and  Force  Components 
From  the  Experimental  Data. 

This  program  determines  the  pressure  distribution  and  calculates 
the  radial  and  tangential  force  components  using  the  experimental 
data.  The  experimental  data  (output  of  the  pressure  transducers 
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In  mv.)  are  input  into  the  program.  Using  the  calibration  equations 

of  the  transducers  (see  Appendix  III)  the  pressure  distribution 

about  the  journal  is  determined.  The  pressure  distribution  is 

then  output  in  ten  degree  increments  of  6 . Using  the  pressui 

data  the  integrals  for  the  force  components  (given  below)  are 

numerically  integrated  using  Simpson's  Rule: 

n2tt 


Fx  = -Lr 


Fy  = -Lr 


■/ 


p sin0  dO 


p cos0  d0 


The  radial  (Fy)  and  tangential  (Fx)  force  components  are  then 
output  from  the  program. 


4.  Calculation  of  the  Pertinent  Squeeze  Film  Parameters  for  the 
Offset  Orbit  Case. 

This  program  vises  the  geometric  relationships  of  the  test  rig  to 
calculate  the  values  of  n , n , 0 , and  a)c  as  functions  of  the 
angular  position  of  the  offset  eccentric  shaft.  The  input  is  the 
dimensions  of  the  test  rig  and  the  thickness  o the  shims  used 
to  offset  the  eccentric  shaft.  The  output  is  n,  n , u>c  , and 
the  position  of  the  line  of  centers,  as  a function  of  the  eccentric 
shaft  position,  given  in  ten  degree  increments  of  the  shaft  rotation. 
The  program  then  uses  the  "long  bearing"  theory  to  predict  the 
pressures  experienced  by  the  pressure  transducer  as  the  eccentric 
shaft  rotates.  The  output  is  the  predicted  pressure  distribution 
as  experienced  by  the  pressure  transducer  in  ten  degree  increments 
of  the  eccentric  shaft  position. 
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5.  Calculation  of  the  Four  Force  Coefficients  for  the  Offset  Orbit 
Case. 

This  program  determines  the  four  force  coefficients,  as  defined 
by  Tondl  [9],  from  the  experimental  data  for  two  positions  of 
journal.  Af  proposed  by  Tondl,  these  force  coefficients  are 
functions  of  n only.  For  the  offset  orbit  case,  there  are  two 
positions  of  the  eccentric  shaft  which  result  in  identical  values 
of  n . Thus  the  force  coefficients  for  these  two  positions 
should  be  identical. 


Fx; 

Fl(nj 

Fy, 

- « M 

F2(n) 

1 

Fx2 

F3(n) 

Jy2 

F4(n) 

3 2 

c = 6tt  y L Tj/c 


0 0 <'2lVi 

(-2wc)1  (2n)x  0 

0 0 (-2 u ) 

C 2. 

(-2u>c)2  (2n)2  0 


(2n)1 

0 

<2n)2 

0 


The  subscripts  ■'  ) .J  (2)  refer  to  the  two  positions  of  the 
eccentric  shaft  with  identical  n . The  force  coefficients  may 
then  be  solved  for: 


Fl(n)” 

■ i [*]” 

F2(n) 

Fyi 

F3(n) 

F»2 

F4(n) 

b i 
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I 


r 
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The  force  components,  Fx^,  Fy^,  Fx2  and  Fy2  are  determined  by 
numerically  integrating  the  pressure  data  for  the  two  journal 
positions.  The  inverse  of  matrix  A is  found  using  the  scientific 
subroutine  MINV.  Matrix  multiplication  is  then  carried  out  to 
deterimine  the  force  coefficients.  The  shaft  position,  position 
of  the  line  of  centers,  n , n ,<*>c  , Fx^  , Fy^  , Fj^  , Fy^  and  the 
force  coefficients  Fl(n)  , F2(n)  , F2(n)  and  F4(n)  are  output. 

/I 


A 


PROGRAM  1 


PRESSURE  O I STR.  IOUT  ION  - • • LONG  BEARING'  • T HEOR  Y 

1 C— • 062 

2 WS=184.94 

3 ED=0. 

4 E = .S 

5 1=0 

6 U=. 0000273d 

7 •’  J=2  • 5 

d .:Pfi=-2. 

9 X=6*U*R J*R J*WS/(C*C  ) 

10  P0=50./X 

11  WRITE  (6,4) 

12  4 FORMAT  (• l './//✓////. T22.* THETA* ,T33, 'PRESSURE • ,T46. 

l'NONJIM,  PRESS.*) 

13  2 IF  (I.GT.360)  GO  TJ  1 

14  TH= 1/67.2  96779 

13  Z=( 14E*C0S(TH>)**2 

16  P=X*(  P04(WPR*E*SIN(  TH)*(24E*C0S(TH)  ) ) / ( ( 2 4E*E  ) *Z  M ( Eu/(  2*h) 
1)*( (l/Z)-(l/(  1 4-E ) * * 2 ) ) *2/W  S) 

17  OPJ=(P-PO*X)*6./X 

Id  TH=TH*S7.295779 

19  WRITE  (6,5)  I .P.OPJ 

20  5 FORMAT  (I25.2F15.4) 

21  1 = 14:10 

22  GO  TO  2 

2T  1 CONTINUE 

24  WRITE  (6,18) 

25  18  FORMAT  ( • 1 • ) 

26  STOP 

27  ENO 


SOATA 


Reproduced  from 
bes)  available  copy. 
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1 

2 

3 

4 

5 

6 
7 
d 
9 

10 

11 

12 


14 

lb 


17 

Id 

19 

20 
21 

22 

23 

24 
2 b 
2d 


PROGRAM  2 


PRESSURE  U iSTRlluUT  ! ON 


••SHORT  BEARING' 


THEORY 


0=. 00002739 
61  - • 5 
E0=0. 

WS=1  B4*<  94 

PU=50. 

EL=4#B13 


C=.062 

1=0 


X3=0. 

RJ=2.5 
WRITE  (6.2) 

2 FORMAT  {• 1 • .✓/////// 
1‘NONDIM.  PPESS.') 

I IF  (I.GT.360)  GO  TO 
TH=I/57. 295779 
Z=1 +&4C0S ( TH ) 
P=P0*6*U* (£0*  COS(  TH) 


,T22 , 'THETA' . T3 3. • PRESSURE* .T46. 
3 

*WS*L*SIN<  TH)  )*(  X3*X3-CL*EL/4*  '( 


C*C4( 


lOP=(P-P0)<‘C*C/(0*WS»RJ*KJ) 
WRITE  (6. 4)1. P. DP 
4 FORMAT  (I25.2F15.4) 

1=1*10 
GO  TO  1 
3 CONTINUE 

WRITE  (6. Id) 

Id  FORMAT  ( ' 1 ' ) 

STOP 

END 


Reproduced  from 
^-st  available  copy. 


2* 


40ATA 


PROGRAM  3 


PRESSURE  Oisr.  AND  FORCE  COMPONENT S - EXPERIMENTAL  OATA 


l 

J l MENS  I ON  P(  J 7)  ,PX(  37)  ,PY(37  ) 

2 

OAT  A P/  8 « * 7 *64*  7. 3,6,92, 6.4, 5. 82, 
l • 94, 2»  82 ,2,4, 2. 6, 3, 34 , 4 t 76  *6*7  8(8 
16.12.4,11.80, 11.3,10.84.10.54,10. 

P0=50. 

4 

RJ  = 2.5 

l. 

U=. 00002 738 

6 

W 5=164. 94 

7 

C=.0b2 

6 

K=0 

9 

H=10/57. 295779 

10 

TH  = 0. 

11 

PTTX=0 

12 

PTTY=0 

1 3 

1 = 1 

14 

J = 1 

15 

WRITE.  (6,12) 

16 

12 

FORMAT  ( • 1*  ./////✓//. T22,*  THETA' . 
1 • NGND 1 M « PRESS.') 

17 

9 

IF  (J.GT.37)  GO  TO  10 

Id 

P( J»=8.l4b*P{ J)-l 1.349 

19 

DP=(P( J )-PU)*C*C/(U*WS*RJ*RJ) 

20 

48 I TE  (6.11)  K.P(J).OP 

21 

1 1 

FORMAT  (I25.2F15.4) 

22 

K=K+10 

23 

3= J + l 

24 

GO  TO  9 

25 

10 

CONTINUE 

2b 

8 

IF  ( I «GT .37)  GO  TO  7 

27 

PXC  l)=P(I )*SIN(TH) 

28 

PY(  I)=-P(I  )*COS(TH) 

29 

1 = 14-1 

30 

TH=TH+H 

31 

GO  TO  8 

32 

7 

CONTINUE 

33 

1 = 1 

34 

c 

IF  (1.GT.36)  GJ  TO  1 

35 

PTX=(PX( l )+4*PX( 1 41  )+PX( 1 + 2) )*H/3 

36 

PTY=(PY( 1 )+4*PY(I +1 ) +PY ( 1+2) )*H/3 

37 

PTTX=PT  TX+PTX 

36 

PTTY=PTTY+PTY 

39 

I = I 4.2 

40 

GO  TO  2 

41 

1 

CONTINUE 

42 

FX=-4.813*2.S*PTTX 

43 

FY=-4.dl3*2.5*PTTY 

44 

WRITE  (6.6) 

45 

6 

FORMAT  ( // X « T 32 , • F ( X ) «.T47.'F(Y)' 

46 

WRITE  (6,5)  FX.FY 

47 

5 

FORMAT  I IF 36. 4,1 F 15. 4) 

4 6 

WRITE  (6,18) 

49 

18 

FORMAT  (•!•) 

50 

STOP 

51 

END 

AO  A TA 


PROGRAM  4 


SQUEEZE  FILM  PARAMETERS  FOR  OFFSET  ORBIT  CASE 


1 

2 

J 

4 

5 
C» 

7 

8 
9 

10 

11 

12 
I J 
14 
lb 


16 

17 

18 

19 

20 
21 
22 

23 

24 
2‘j 
26 

2 7 
23 

29 

30 

31 

32 
23 

34 

35 

36 

3 7 

38 

39 

40 

41 

42 

43 

44 
4b 

46 

47 
4M 

49 

50 
61 
u d 

»J 

64 

'»■> 

‘>6 

b/ 


68 

59 

60 
•>  l 
62 
6 ? 
04 

6 b 
60 
6 7 

<>M 
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DIMENSION  ENI 37 ) « fcNO 137) * WCC( 371 , THE (37) 

I  =0 

S=.0I5 
*S=184.94 
E P= • 031 
GA=0.0 
C=.062 
J = 1 

THO=20. 

U=,0000273u 

RJ=2.5 

X=6*U4R J*R3*^S/(C»C» 

P(J=40/X 
WRITE  (6,2) 

2 FORMAT  (• 1* .✓/✓// ///.T18,» SHAFT  ANGLE  * ,T  36  •»  N*  , T46,  • N-DOT  ♦ , 
IT 52.* THETA* ,T74,«  W-C* ) 

3 IF  ( I .GT.180)  GO  TO  5 

6  — SQRT(  S*S*2*S*eP  *CQS( GA ) fEPOEP) 

EO=-( S*EP*SIN(GA> *W S/E) 

IF  (I .EQ.O)  GO  TO  6 

IF  (I.LQ.180)  GO  TO  7 

Z=( (S*S)F(t*E)-(£P*EP) )/( 2 ,*S*E ) 

TH=ARCOS(Z) 

wC=-(  (ED/S)-(Z*EO/E)l/SGRTU.-(Z*Z)  ) 

GO  TO  9 
o TH=0.0 
GO  TO  8 

7 TH=180/57.2958 

8 wC=EP*wS/fe 

9 CONTINUE 
14  = 38“  J 
EN(J)=E/C 
EN(  K) =EN( J ) 

ENO ( J ) = ED/C 
ENO(K) =-END( J ) 

WCC( J ) =*C 
WCC(K  ) =WC 

THE( J)  = TH*S7.29b7  79 
THE(K) =360,— THE (31 
J=J*1 
1 = 1*10 

GA  = I/b7. 295779 
GO  TO  3 
b CUNTINUE 
GA=0.0 
J=1 

10  IF  (J.GT.37)  GO  TO  11 

WRITE  (6,12)  GA»E  N( J ) » END (J)  » THE { J ) , WCC( J ) 

12  FORMAT  ( 1F2S.  1 . 2F  1 4 .3,  iF*  12 . 1 .1F13.2) 

J = JF1 
GA  = GA«-10 
GO  TO  10 

11  CONTINUE 
G A = 0 • 0 

I =0 
J = 1 

WRITE  (6,13) 

12  FORMAT  C • 1 •,////////, r\0 ,* SHAFT  ANGLE • »T 36, * THETA *, T50 , 

1* PRESSURE* ) 

14  IF  (i.GT.JOO)  GO  TO  lb 
TH=THUM80-THc(  J) 

IF  ( TH,GT  ,0  ) GO  TC  1j 
TH=3oO*  TH 
16  CONTINUE 

TH  = TH/5  7«  29S  7 79 
E=EN(  j) 

f:no(  j)=enu(  j ) *2/wcct  j) 

WPR  = - (2**CL(J)/WS > 

Z=( IFEN( J)*COS( TH) )**2 

P = X*  ( POF (WPR*E*SIN( TH)* (2FEN ( j ) *COS (TH) ) )/( (2*EN(J)*EN( J) I* 
1 Z ) + ( END  ( J )/  ( 2 *EN(  0)  ))*(  (l/Z)“(  l/(  ’.FEN(  J ) ) 4*2  ) I * ('2^WS)  ) 
TH=TH*57. 295779 
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o ~ ''J  C *W  ?>  > O 


f } *\U  Or  Am 

r i.<(,r  Cut_>  r-  ic  i * nt  ii  - l k>~ o-  imcnt  al.  oata,  oh‘:>(  t u-i>it  casc 

i>  I *•!*.  Nit  ON  A ( A . 4 ) . ,.N  o ( ) , Fi*  ( 4 > , F X ( ; ) , F Y ( Z ) , L ( . ) , M { 4 ) 

C - •Oo>”; 

c:  •_  =4 . a 1 3 

i-.OlD 
•V  i = I U4  .04 
*:  f’  -.OJ  1 
U=  *JOO 02730 

X-u‘U«r.L*t.U*>*3)/(C‘C) 

C»A=  <0/  a?»2  -Id? /y 
:-N=r./C 

cl)  = -(D*tP*S[N(  GA)  »*#  S ) /t: 

CNL)(  1 > =tD/C 
LND(  2)  =-£ND(  1 ) 

^ = ( :.<'S+t.*d-hP*t:P)/(  2*3*fcJ 
PH  I -=Ai-lC03C  Z ) 

*C  =-(£U/S~Z»ED/£)/SQ«T (1-Z*Z> 

(6*15) 

1 bi T61?* 5-c • ) * ,//////// ’Tc0 • r 31 . • N-DOV ( 1 ) • . r^o.  • N-0UT(  Z) • . 

(o,16>  LN.CNiHl  ),UND(2)  ,*C 
16  FORMAT  ( 1F22.J, 3r 15.3) 

A<  1 .1  )=G. 

A(  1 ,2) =0. 

A { 1 .3) =-2**C 
A(  1,4)  =2*EN0(  1 ) 

A(  2.1  ) --?_** C 
A( Z. 2)  -2-*CNi>{  1 ) 

A(  2.  j)  =u. 

A(  .2*41=0* 

A(  3 , l ) =0  • 

A(  J.2) -0. 

A(  3, J) =— 2*  »C 
A(  3*  4 ) =2  *«£NN(  z ) 

A ( 4 , 1 ) =-2*WC 
A<  4*2  > =2  *EN!)(  2 ) 

A{ 4,3) —0  • 

A<  4 ,4) =0. 

N-s4 

CAUL  FORCES  ( F X , f Y ) 

CAUL  M INV(  A»N  *i>  ,L  »M) 

J=1 

li  ^F  (J.GT.4)  oU  TO  10 

J=J>l  = (A<J*l>  FX<  l,*A( J»2)*FY( l )fA( j,  "FX(Z  )FA(J  *4)  *rV(L*  ) )/X 
GO  TO  11 

10  continue 

W-.  ITL  (0.13) 

1 1 *°*  ,//V*T-0  .’FN(  1)*  ,rj>,»FN(Z)  • ,ToO.*FN(3)*  *T63. 

1 7 r N|  4 J 1 J 

V.<ITi_  (0.12)  FN(  1 ) » F N ( Z ) *FN(3)  *(-N(4) 

1 e Pi.cMA)  (IF 24  * 3 * 3F  13*3) 

I It.  (o*20) 

2 0 r-0  r..A  AT  ( • 1 • ) 


i M.i-’-r-j"  ( r X * K Y ) 

’ ■ ir  * ("■.  >*x  (,•),  fy  ( 2 ),'jar)  ,px(2  n 

~ f-  r/i  ,'?.i  7,2.  36, 2. 

« 1 • * JS  : 1 • jj  « • •#V.4  * • 6?  i • Y7  , X 


t V * , 1 . < « : 1 . 3 
! l , 

.HM  0/1  .9  3 

i /- 1 . . i . r.  > 

l :• . ? , ;* . $-5,  .> 

K-  i 

»'•-*  7.2*?*>7V9 


, • • #V.4  « • 6?  i • Y7  , x 

• tv*’ » • 'J » • y i>  ,/*  • 1 • 

i .o».i  .7?.  i .o,  i . :*,  t , 

v.  22, ji  , 3. iiii  ,3. 3* . j. 
•f-i.  2.22.  2.2  7,2.2-3.2.27 


t L'Y( 
30,2. 

m * * * < 

l.  / /, 
. :lr  , , 
3 7.3, 
* 2.2. 


. o ( 2 r ) 

> * ?-  • 1 .5  • 
,•*.  »r;,  i 
. n i,  i . « 

i ■■>.- 
,.(.07  ,2 
>.  1 7,2. 


1 . O'-  , 1 
.0  «. 
••1.1  . >*• 


! t r~ 

•I'lftr" 

J = 

1 

t *; 

( r - ^ ) 0.1 

<y  r "s 

*.*■  I MU' 

T t* 

( J.«T.  37) 

:>(  J >=-.*«.  1'--  *-M  J 
J-  Jfl 
or  to  c, 

0 '■'**■'  IMU'* 

«r  (J.*'.T.37)  05  T 
1 c J)  -g.*m. ; *: o ( j 

j-  j ■*■  l 
0?  - •»  tft 

’ our 

d?(). 

! =1 

ik  (r.iT.?7)  nn  t 
»X  M )-P(  ! )*$!N(  G) 
r»Y(  : ) -p(  r ) *C03(G) 
t = *f  l 
0=0  *-H 

00  ro  ?> 

C’N”  1 7|IJT 


•;o  ro  13 

T->CJ)  — 11.  3^9 


G‘i  TO  1 1 

•=o(  j)-i  i . :uo 


T = 1 
PT  t 

pt  ’ 

T 

T .< 

,>-Y 

L*T  ~ 
r>T-: 

i= : 
c,~ 

C 

r X ( 
rY  ( 

V<‘-  ! 
i;  -i 
,»•  7 

j ... 

v r) 


X =0  . 

Y=0. 

( r • G‘r.?'  ) O'*  TO  ! 

-<'X<  l )*'  *.»X<  I 4-1  ) f>X(  1*2)  >*H/3. 
= |1'Y(  : 1 *■««!' Y(  [ 1 ) >PY(  ! ft*)  ) »M/3. 
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